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Nanofluids  are  considered  to  have  great  potential  for  heat  transfer  enhancement  and  are  highly  suited 
to  application  in  practical  heat  transfer  processes.  Recently,  several  important  studies  were  carried  out 
to  understand  and  explain  the  causes  of  the  enhancement  or  control  of  heat  transfer  using  nanofluids. 
The  main  aim  upon  which  the  present  work  is  based  is  to  give  a  comprehensive  review  on  the  research 
progress  on  the  natural  convective  heat  transfer  characteristics  of  nanofluids  for  both  single-  and  two- 
phase  models.  Both  experimental  and  theoretical  studies  are  reviewed  for  natural  convection  of 
nanofluids  in  different  types  of  enclosures. 
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1.  Introduction 

Natural  convection  heat  transfer  is  an  important  phenomenon  in 
engineering  and  industry  with  widespread  application  in  diverse 
field,  such  as,  geophysics,  solar  energy,  electronic  cooling,  nuclear 
energy  and  other  industrial  system  in  various  sectors.  A  major 
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limitation  against  enhancing  the  heat  transfer  in  such  engineering 
systems  is  the  inherently  poor  thermal  conductivity  of  conventional 
fluids,  including  oil,  water  and  ethylene  glycol  mixture.  Therefore,  for 
more  than  a  century  since  Maxwell’s  theories  in  1873,  scientists  and 
engineers  have  made  great  effort  to  break  this  fundamental  limit  by 
dispersing  millimeter-or  micrometer  sized  particles  in  liquids.  How¬ 
ever,  the  major  problem  with  the  use  of  such  large  particles  is  the 
rapid  sedimentation  of  these  particles  in  fluids. 

Maxwell’s  concept  is  old  but  what  is  new  and  innovative  in  the 
concept  of  nanofluids  is  the  idea  of  using  nanometer-sized 
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Nomenclature 

T 

stress  tensor  (Pa) 

a 

thermal  diffusivity  (m2  s-1) 

cp 

specific  heat  at  constant  pressure  (J  kg-1  K-1) 

</> 

dimensional  nanoparticle  volumetric  fraction 

Dj 

thermal  diffusion  coefficient  (m2/s) 

</> 

inclination  angle 

d 

diameter  (m) 

p 

density  (kg  m-3) 

g 

gravitational  acceleration  (ms-2) 

n 

ratio  of  the  nanolayer  thickness  to  the  original 

H 

dimensional  cavity  height  (m) 

particle  radius 

k 

thermal  conductivity  (W/m  K) 

1<B 

Boltzmann  constant  (J/K) 

Subscripts 

P 

pressure  (Pa) 

Pr 

Prandtl  number,  Pr  =  Vf/oCf 

/ 

base  fluid 

Pe 

Peclet  number 

nf 

nanofluid 

Re 

Reynolds  number,  Re  =  pfkBT/3TC/ijlf 

P 

particle 

T 

dimensional  temperature  (°C) 

k 

thermal  conductivity  (W/m  K) 

V 

dimensional  nanofluid  velocity  (ms-1) 

Greek  symbols 

P 

dynamic  viscosity  (Ns  m  2) 

particles  to  create  stable  and  highly  conductive  suspensions, 
primarily  for  suspension  stability  and  for  dynamic  thermal  inter¬ 
actions.  Recognizing  an  excellent  opportunity  to  apply  nanotech¬ 
nology  to  thermal  engineering,  Eastman  et  al.  [1]  is  the  first  who 
conceived  the  novel  concept  of  nanofluid  by  hypothesizing  that  it 
is  possible  to  break  down  these  century-old  technical  barriers  by 
exploiting  the  unique  properties  of  nanoparticles. 

Nanofluids  are  dilute  suspensions  of  nanoparticles  smaller  than 
100  nm.  The  thermal  conductivity  of  the  particle  materials,  metal¬ 
lic  or  nonmetallic  such  as  A1203,  CuO,  Cu,  SiO  and  Ti02  are 
typically  order-of-magnitude  higher  than  the  conventional  fluid. 
Such  nanofluids  represent  a  class  of  its  own  different  from  those  of 
traditional  solid-mixtures,  they  have  superior  properties  like  high 
thermal  conductivity  [1-5],  minimal  clogging  in  flow  passage, 
reduced  puming  power,  long  term  stability  and  homogeneity. 
Therefore  nanofluids  have  a  wide  range  of  potential  application 
where  improved  heat  transfer  is  required. 

The  purpose  of  this  paper  first  is  to  review  the  available  results  of 
the  flow  and  heat  transfer  presented  in  literature  for  natural 
convection  of  nanofluids,  including  the  experimental  and  theoretical 
investigations  for  both  single-  and  two-phase  models.  The  results  of 
massive  number  of  different  parameters  affect  the  fluid  flow  and 
heat  transfer  characteristics  are  summarized.  The  second  purpose  is 
to  explore  the  causes  of  the  enhancement  of  heat  transfer,  which 
would  be  of  help  in  understanding  the  transport  in  nanofluids,  and 
accurately  predict  natural  convection  using  nanofluids. 


2.  Brief  theory 

2.2.  Model  description 

In  the  literature,  convective  heat  transfer  with  nanofluids  can 
be  modeled  using  mainly  the  two-phase  or  single  approach.  In  the 
two-phase  approach,  the  velocity  between  the  fluid  and  particles 
might  not  be  zero  [6]  due  to  several  factors  such  as  gravity, 
friction  between  the  fluid  and  solid  particles,  Brownian  forces, 
Brownian  diffusion,  sedimentation  and  dispersion.  In  the  second 
approach,  the  nanoparticles  can  be  easily  fluidized  and  therefore, 
one  may  assume  that  the  motion  slip  between  the  phases,  if  any 
would  be  considered  negligible  [7].  The  latter  approach  is  simpler 
and  more  computationally  efficient. 


2.2.  Conservation  equations 
2 2.1.  Single  phase  model 

Although  nanofluids  are  solid-liquid  mixtures,  the  approach 
conventionally  used  in  most  studies  of  natural  convection  handles 
the  nanofluid  as  a  single-phase  (homogenous)  fluid.  In  fact,  due  to 
the  extreme  size  and  low  concentration  of  the  suspended  nanopar¬ 
ticles,  the  particles  are  assumed  to  move  with  same  velocity  as  the 
fluid.  Also,  by  considering  the  local  thermal  equilibrium,  the  solid 
particle-liquid  mixture  may  then  be  approximately  considered  to 
behave  as  a  conventional  single-phase  fluid  with  properties  that  are 
to  be  evaluated  as  functions  of  those  of  the  constituents. 

The  governing  equations  for  a  homogenous  analysis  of  natural 
convection  are  continuity,  momentum,  and  energy  equations 
with  their  density,  specific  heat,  thermal  conductivity,  and  visc¬ 
osity  modified  for  nanofluid  application.  The  specific  governing 
equations  for  various  studied  enclosures  are  not  shown  here  and 
they  can  be  found  in  different  references  [18,32,33,55]. 

2.2.2.  Two-phase  model 

Several  authors  have  tried  to  establish  convective  transport 
models  for  nanofluids.  Nanofluid  is  a  two-phase  mixture  in  which 
the  solid  phase  consists  of  nano-sized  particles.  In  view  of  the 
nanoscale  size  of  the  particles,  it  may  be  questionable  whether 
the  theory  of  conventional  two-phase  flow  can  be  applied  in 
describing  the  flow  characteristics  of  nanofluid.  On  the  other  hand, 
several  factors  such  as  gravity,  friction  between  the  fluid  and  solid 
particles  and  Brownian  forces,  the  phenomena  of  Brownian  diffusion, 
sedimentation,  and  dispersion  may  affect  a  nanofluid  flow.  Conse¬ 
quently,  the  slip  velocity  between  the  fluid  and  particles  cannot  be 
neglected  for  simulating  nanofluid  flows.  Since  the  two  phase 
approach  considers  the  movement  between  the  solid  and  fluid 
molecule,  it  may  have  better  prediction  in  nanofluid  study.  To  fully 
describe  and  predict  the  flow  and  behavior  of  complex  flows, 
different  multiphase  theories  have  been  proposed  and  used.  The 
large  number  of  published  articles  concerning  multiphase  flows 
typically  employed  the  Mixture  Theory  to  predict  the  behavior  of 
nanofluids  [8-10]. 

A  comprehensive  survey  of  convective  transport  in  nanofluids 
was  made  by  Buongiorno  [11],  using  a  model  in  which  Brownian 
motion  and  thermophoresis  are  accounted  for.  Buongiorno 
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developed  a  two-component  four-equation  non-homogeneous 
equilibrium  model  for  mass,  momentum,  and  heat  transfer  in 
nanofluids.  The  nanofluid  is  treated  as  a  two-component  mixture 
(base  fluid  +  nanoparticles)  with  the  following  assumptions: 

•  No  chemical  reactions. 

•  Negligible  external  forces. 

•  Dilute  mixture  (0  =  1). 

•  Negligible  viscous  dissipation. 

•  Negligible  radiative  heat  transfer. 

•  Nanoparticle  and  base  fluid  locally  in  thermal  equilibrium. 

Invoking  the  above  assumptions,  the  following  equations 
represent  the  mathematical  formulation  of  the  non-  homogenous 
single  phase  model  for  the  governing  equations  as  formulated  by 
Buongiorno  [11]: 

2.2.2. 2.  Continuity  equation. 

Vv  =  0  (1) 


22.2.2.  Nanoparticle  continuity  equation. 

^v-v0  =  v-  (dbv4>+dt'22j  (2) 


Here  0  is  nanoparticle  volume  fraction,  DB  is  the  Brownian 
diffusion  coefficient  given  by  the  Einstein-Stokes’s  equation: 


where  p  is  the  viscosity  of  the  fluid  and  dp  is  the  nanoparticle 
diameter  and  DT  is  the  thermophoretic  diffusion  coefficient, 
which  is  defined  as 


k 

2  k+kp 


2.3.  Physical  properties  of  the  nanofluids  for  single-phase  model 

Base  nanofluid  properties  have  been  published  over  the  past 
few  years  in  literature.  However,  only  recently  have  some  data  on 
temperature-dependent  properties  been  provided,  even  though 
they  are  only  for  nanofluid  effective  thermal  conductivity  and 
effective  absolute  viscosity. 

2.3. 2.  Density 

In  the  absence  of  experimental  data  for  nanofluid  densities, 
constant-value  temperature  independent  values,  based  on  nano¬ 
particle  volume  fraction,  are  used: 

Pnf  =  0  —  0)P/  +  0 Pp 


2.3.2.  Specific  heat  capacity 

It  has  been  suggested  that  the  effective  specific  heat  can  be 
calculated  using  the  following  equation  as  reported  in  Refs. 
[12-15]  as 

(Cp)n/  =  (1  ~ 0)(Cp)/  +  0(Cp)p 

Other  authors  suggest  an  alternative  approach  based  on  heat 
capacity  concept  [7,16]: 

(pCp)nf  =  (1  ~  0)(pCp)/  +  (t>(pCp)p 

These  two  formulations  may  of  course  lead  to  different  results 
for  specific  heat.  Due  to  the  lack  of  experimental  data,  both 
formulations  are  considered  equivalent  in  estimating  nanofluid 
specific  heat  capacity  [14]. 

2.3.3.  Dynamic  viscosity  and  thermal  conductivity 

In  this  work,  different  nanofluid  models  based  on  a  combina¬ 
tion  of  the  different  formulas  for  the  dynamic  viscosity  and 
thermal  conductivity  adopted  in  the  studies  of  natural  convection 
are  summarized  in  Tables  1  and  2. 


In  Eq.  (5),  k  and  kp  are  the  thermal  conductivity  of  the  fluid  and 
particle  materials,  respectively. 


3.  Numerical  investigations 


2.22.3.  Momentum  equation. 

v-Vv  = — —  Vp  +  V-x+g  (5) 

Pnf 

where 

T=-/'n/(VF  +  ^VF)t)  (6> 

where  the  superscript  ‘t’  indicates  the  transpose  of  Vv. 

2.22.4.  Energy  equation. 

v-VT=  V(ocnfVT)  +  E22-  (dbV4>  VT +Dj  SLYTj  (7) 

This  nanofluid  model  can  be  characterized  as  a  ‘two-fluid’ 
(nanoparticles  +  base  fluid),  four-equation  (2  mass  +  1  momentum+ 
1  energy),  non-homogeneous  (nanoparticle/fluid  slip  velocity 
allowed)  equilibrium  (nanoparticle/fluid  temperature  differences  not 
allowed)  model. 

Note  that  the  conservation  equations  are  strongly  coupled. 
That  is,  v  depends  on  0  via  viscosity;  0  depends  on  T  mostly 
because  of  thermophoresis;  T  depends  on  0  via  thermal  con¬ 
ductivity  and  also  via  the  Brownian  and  thermophoretic  terms  in 
the  energy  equation:  0  and  T  obviously  depends  on  v  because  of 
the  convection  terms  in  the  nanoparticle  continuity  and  energy 
equations,  respectively. 


3.2.  Single  phase  model 

Khanafer  et  al.  [17]  were  the  first  who  studied  numerically,  the 
buoyancy  driven  heat  transfer  enhancement  of  nanofluid  in  a 
square  cavity.  They  used  Wasp’s  model  [18]  for  thermal  con¬ 
ductivity  and  Brinkman  model  [19]  for  viscosity.  They  illustrated 
that  the  suspended  nanoparticles  substantially  increase  the  heat 
transfer  rate  at  any  given  Grashof  number.  In  addition,  they 
indicated  that  the  nanofluid  heat  transfer  rate  increases  with  an 
increase  in  the  nanoparticles  volume  fraction.  They  also  showed 
that  the  presence  of  nanoparticles  in  the  fluid  is  found  to  alter  the 
structure  of  the  fluid  flow. 

Jang  and  Choi  [20]  examined  the  thermal  characteristics  of  free 
convection  in  rectangular  cavity  heated  from  below  with  nano¬ 
fluid  such  as  water  based  nanofluids  containing  6  nm  copper  and 
2  nm  diamond,  using  Jang  and  Choi’s  model  [21]  for  the  effective 
thermal  conductivity.  They  concluded  that  the  nanofluids  is  more 
stable  and  have  larger  heat  transfer  coefficient  compared  with  a 
classical  fluid  in  a  rectangular  cavity  with  free  convection. 

Kim  et  al.  [22]  analyzed  the  convective  instability  driven  by 
buoyancy  and  heat  transfer  characteristics  of  nanofluids.  They 
chose  Wang  et  al.  [23]  model  based  on  the  mean  field  approach 
for  expressing  the  thermal  conductivity  enhancement.  They 
indicated  that  as  the  density  and  heat  capacity  of  nanoparticles 
increase  and  the  thermal  conductivity  and  the  shape  factor  of 
nanoparticles  decrease,  the  convective  motion  in  a  nanofluid  sets 
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Table  1 

Various  viscosity  models  for  nanofluids. 

Models 

Brinkman  model  [19] 

Einstein  model  [27] 

Pak  and  Cho’s  Correlation  [28] 

Nguyen  et  al.  model  [36]  (fitted  by  Abu-Nada  [34]) 


Jang  et  al.  model  [41] 


Koo  and  Kleinstreuer  [48] 


Maiga  model  [81  ] 

Brownian  model  [82] 
Nguyen  model  [83] 
Gherasim  et  al.  [84] 


Dynamic  viscosity 

[if 

[inf  =  /*/(l  +  2.5(f))  (f)  +  0.05 
[ Lnf  =  [if(  1  +39.1 1(/>  +  533.9(/>2) 


ncnc_  15.937  .  ,  1356.14 

Hciio  —  — 0.6967  -f - — - (- 1 .238  (f)  -f - ^ - 


-0.259(/)2-30.88^-196^3274+0.01593(/)3 +4.38206^ +  147.573^ 


[iAl  0  =  exp(3.003—0.04203T— 0. 54450 +  0.0002553T2  +  0.0524(/)z  — 1.622  0-1) 


[inf  —  [if  ('l  +2.5  cf)) 


—2s 


1  (tt)  ^2/3(£+1) 


flB  =  5x1 J  (ffff  0.4) 

/3  =  0.01 37(1 00(f))-0 8229  for  (f)<  1% 

(3  =  0.001 1(1  OO^r0 7272  for  (f)>\% 

f(T,(f))  =  (-6.04(f)  +  0.4705)1+ (1 722.30-134.63) 
1%  <  (f)  <  4%  and  300  k  <  T  <  325  k 

[inf  =  [if(f  +  7.3(f)  + 123<^>  ) 

iunf  =  A/(l  +  2..5(f)  +  6.1 70  ) 

[inf  =  /^(l  +0.025(/>  +  0.015(/)2) 

=  0.904eU8^ 


in  easily.  They  also  showed  that  heat  transfer  coefficient  of  a 
nanofluid  is  enhanced  by  all  parameters  with  respect  to  the 
volume  fraction. 

Jou  and  Tzeng  [24]  numerically  studied  the  natural  convection 
heat  transfer  enhancements  of  nanofluids  within  a  two-dimensional 
enclosure.  They  analyzed  heat  transfer  performance  using  Khana- 
fer’s  model  for  various  parameters,  like  volume  fraction,  Grashof 
number,  and  aspect  ratio  of  the  enclosure.  Results  showed  that 
increasing  the  buoyancy  parameter  and  volume  fraction  of  nano- 
fluids  cause  an  increase  in  the  average  heat  transfer  coefficient. 

Polidori  et  al.  [25]  analytically  investigated  the  natural  con¬ 
vection  heat  transfer  of  Newtonian  nanofluids  over  a  vertical 
semi-infinite  plate.  In  their  study,  the  dealt  with  y-Al203/water 
nanofluids  whose  Newtonian  behavior  was  experimentally  con¬ 
firmed  for  particle  volume  fractions  less  than  4%.  It  was  shown 
that  special  care  has  to  be  exercised  in  drawing  generalized 
conclusions  about  the  heat  transfer  enhancement.  They  found 
that  natural  convection  heat  transfer  is  not  solely  characterized 
by  the  nanofluid  effective  thermal  conductivity  and  that  the 
sensitivity  to  the  viscosity  model  used  seems  undeniable  and 
play  a  key  role  in  the  heat  transfer  behavior. 

Hwang  et  al.  [26]  theorically  investigated  the  thermal  char¬ 
acteristics  of  natural  convection  in  a  rectangular  cavity  heated 
from  below  with  water-based  nanofluids  containing  alumina 
(A1203).  They  used  thermal  conductivity  model  of  Jang  and  Choi 
[21]  and  various  models  for  the  effective  viscosity.  It  was  shown 
that  the  experimental  results  are  put  between  two  theoretical 
lines  derived  from  Jang  and  Choi’s  conductivity  model  used  along 
with  both  Einstein’s  model  [27]  and  Pak  and  Cho’s  correlation  [28] 
for  viscosity.  They  showed  that  water-based  A1203  nanofluids  is 


more  stable  than  base  fluid  in  a  rectangular  cavity  heated  from 
below  as  the  volume  fraction  of  nanoparticles  increases,  the  size 
of  nanoparticles  decreases,  or  the  average  temperature  of  nano¬ 
fluid  increases.  They  also  showed  that  the  ratio  of  heat  transfer 
coefficient  of  nanofluids  to  that  of  base  fluid  is  decreased  as  the 
size  of  nanoparticles  increases,  or  the  average  temperature  of 
nanofluids  is  decreased.  In  addition,  they  stated  that  the  effective 
thermal  conductivity  models  used  by  Kim  et  al.  [22]  and  Khanafer 
[17]  cannot  predict  the  effect  of  the  size  of  nanoparticles  and  the 
average  temperature  of  nanofluids  on  the  enhancement  of  the 
effective  thermal  conductivity,  although  the  models  are  able  to 
predict  the  effect  of  the  particle  fraction  of  nanoparticles  on  the 
enhancement  of  the  effective  thermal  conductivity  of  nanofluids. 

Ho  et  al.  [29]  aimed  to  clarify  and  present  a  deep  explanation 
on  the  influences  of  uncertainties  due  to  adopting  various 
formulas  for  the  effective  thermal  conductivity  and  dynamic 
viscosity.  They  presented  a  numerical  analysis  for  natural  con¬ 
vection  of  water  A1203  nanofluid  in  a  square  vertical  enclosure. 
They  observed  that  a  significant  difference  between  enhance¬ 
ments  in  the  dynamic  viscosity  estimated  from  the  two  adopted 
formulas  leads  to  contradictory  heat  transfer  efficacy  of  the 
nanofluid,  so  that  the  heat  transfer  across  the  enclosure  can  be 
found  to  be  enhanced  or  migtated  with  respect  to  the  base  fluid. 
According  to  the  authors,  the  contradictory  behavior  dictated  by 
the  significant  difference  between  the  dynamic  viscosity  enhance¬ 
ments  calculated,  using  two  formulas  might  contribute  to 
explaining  the  disparate  findings  among  the  numerical  predic¬ 
tions  [17]  and  experimental  results  [76,77]. 

Oztop  and  Abu-Nada  [30]  investigated  heat  transfer  and  fluid  flow 
due  to  buoyancy  forces  in  a  partially  heated  enclosure  using  different 
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Table  2 

Various  thermal  conductivity  models  for  nanofluids. 


Model 


Thermal  conductivity 


Wasp  model  [18] 
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Jang  and  Choi  model  [21] 
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Bruggeman  model  [23] 


Chon  et  al.  model  [35] 
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1  +64.70°  7640  pr0.9955Re1.2321 


Koo  and  Kleinstreuer  [48] 


Charuyakorn  et  al.  [85] 


Staionary  model  [86] 


Yu  and  Choi  [87] 
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Patel  et  al.  [88] 


Mintsa  et  al.  [89] 
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nanofluids.  The  flush  mounted  heater  was  located  to  the  left  vertical 
wall  with  a  finite  length.  The  temperature  of  vertical  wall  was  lower 
than  that  of  heater  while  other  walls  were  insulated.  Calculations 
were  performed  for  different  Rayleigh  numbers,  height  of  heater, 
location  of  heater,  aspect  ratio  and  volume  fraction  of  nanoparticles. 
They  found  that  heat  transfer  increases  with  increasing  Rayleigh 
number,  height  of  heater  and  volume  fraction  of  nanoparticles.  In 
addition,  the  heat  transfer  enhancement  was  found  to  be  dependent 
on  the  type  of  nanofluid  and  more  pronounced  at  low  aspect  ratio 
that  at  high  aspect  ratio. 

Ogiit  [31]  studied  numerically  the  natural  convection  heat 
transfer  of  water-based  nanofluids  in  an  inclined  enclosure  with  a 
constant  flux  heater.  Five  types  of  nanoparticles  were  taken  into 
consideration:  Cu,  CuO,  A1203  and  Ti02.  The  results  showed  that 
the  average  heat  transfer  rate  increases  significantly  as  Rayleigh 
number  and  particle  volume  fraction  increase,  the  effect  is  more 
pronounced  for  Ag  and  Cu  nanoparticles.  The  results  also  showed 
that  the  average  heat  transfer  decreases  with  an  increase  in  the 
length  of  the  heater.  As  the  heater  length  is  increased,  the  average 
heat  transfer  rate  starts  to  decrease  for  a  smaller  inclination  angle. 

Abu-Nada  and  Oztop  [32]  analyzed  numerically  the  effect  of 
inclination  angle  on  natural  convection  heat  transfer  and  fluid  flow 
in  a  two  dimensional  enclosure  filled  with  Cu-nanofluid.  It  was  found 
that  the  effect  of  nanoparticle  concentration  on  Nusselt  number  is 
more  pronounced  at  high  volume  fraction  than  at  low  volume 
fraction.  It  was  also  found  that  the  Effects  of  inclination  angle  on 
percentage  of  heat  transfer  enhancement  become  insignificant  at  low 
Rayleigh  number  but  it  decreases  the  enhancement  of  heat  transfer 
with  nanofluid. 

Abouali  and  Falahatpisheh  [33]  carried  out  a  numerical  study 
of  internal  free  convection  of  A1203  water  nanofluid  in  vertical 
annuli.  Vertical  walls  are  maintained  at  constant  temperatures 


and  horizontal  walls  are  adiabatic.  They  developed  a  correlation 
as  a  function  of  Nusselt  number  of  base  fluid  and  particle  fraction 
which  is  a  monotonically  linear  decreasing  function  of  particle 
fraction.  In  addition,  they  showed  that  thermal  conductivity 
model  of  Wasp  [18]  greatly  overestimates  the  Nusselt  number 
of  experimental  data  of  Putra  et  al.  [76]  and  the  thermal 
conductivity  of  Jang  and  Choi  [21]  agrees  well  with  the  experi¬ 
ment  and  shows  deterioration  of  natural  convection  heat  transfer. 

Abu-Nada  [34]  studied  the  heat  transfer  enhancement  in 
horizontal  annuli  using  variable  properties  of  Al203-water  nano¬ 
fluid.  He  used  Chon  et  al.  [35]  expression  for  conductivity  and 
Nguyen  et  al.  [36]  experimental  data  for  viscosity.  They  observed 
that  for  Ra>104,  the  average  Nusselt  number  was  reduced  by 
increasing  the  volume  fraction  of  nanoparticles.  However,  for 
Ra  =  103,  the  average  Nusselt  number  increased  by  increasing  the 
volume  fraction  of  nanoparticles.  For  Ra  =  104,  the  Nusselt  num¬ 
ber  was  deteriorated  everywhere  around  the  cylinder  surface 
especially  at  high  expansion  ratio.  However,  this  reduction  is  only 
restricted  to  certain  regions  around  the  cylinder  surface  at 
Ra  =  103.  They  observed  also  that  The  Nguyen  et  al.  data  and 
Brinkman  model  gives  completely  different  predictions  for 
Ra  =  104  where  the  difference  in  prediction  of  Nusselt  number 
reached  30%.  However,  this  difference  was  less  than  10%  at 
Ra  =  103. 

Aminossadati  and  Ghasemi  [37]  presented  a  study  of  natural 
convection  in  a  partially  heated  enclosure  from  below  and  filled 
with  different  types  of  nanofluids.  The  effects  on  the  enclosure 
cooling  performance  of  Rayleigh  number,  solid  volume  fraction, 
heat  source  length  and  location  and  the  type  of  nanofluid  were 
studied.  The  results  indicated  that  adding  nanoparticles  into  pure 
water  improves  its  cooling  performance  especially  at  low  Ray¬ 
leigh  numbers.  The  type  of  nanoparticles  and  the  length  and 
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location  of  the  heat  source  proved  to  significantly  affect  the  heat 
source  maximum  temperature. 

Abu-Nada  et  al.  [38]  presented  a  study  of  the  heat  transfer 
enhancement  in  a  differentially  heated  enclosure  using  variable 
thermal  conductivity  and  variable  viscosity  of  Al203-water  and 
CuO-water  nanofluids.  The  results  indicated  that  for  an  enclosure 
with  unity  aspect  ratio,  the  average  Nusselt  number  of  A1203- 
water  nanofluid  at  high  Rayleigh  number  was  reduced  by  increas¬ 
ing  the  volume  fraction  of  nanoparticles  above  5%.  However,  at 
low  Rayleigh  number,  the  average  Nusselt  number  was  slightly 
enhanced  by  increasing  the  volume  fraction  of  nanoparticles. 
Also,  at  high  Rayleigh  numbers,  Cuo-water  nanofluids  manifest 
a  continuous  decrease  in  Nusselt  number  as  the  volume  fraction 
of  nanoparticles  is  increased.  However,  the  Nusselt  number  was 
not  sensitive  to  the  volume  fraction  at  low  Rayleigh  numbers.  It 
was  observed  that  enclosures,  having  high  aspect  ratios,  experi¬ 
ence  more  deterioration  in  the  average  Nusselt  number  when 
compared  to  enclosures  having  low  aspect  ratios.  It  was  also 
observed  that  the  average  Nusselt  number  was  more  sensitive  to 
the  viscosity  models  than  to  the  thermal  conductivity  models. 

Lin  and  Violi  [39]  analyzed  the  heat  transfer  and  fluid  flow  of 
natural  convection  in  a  cavity  filled  with  Al203/water  nanofluid 
that  operates  under  differentially  heated  walls.  The  thermal 
conductivity  and  dynamic  viscosity  of  the  nanofluid  are  employed 
by  Xu’s  model  [40]  and  Jang’s  et  al.  model  [41],  respectively.  The 
heat  transfer  rates  were  examined  for  parameter  of  non-uniform 
nanoparticle  size,  mean  nanoparticle  diameter,  nanopartticle 
volume  fraction,  Prandtl  number  and  Grashof  number.  They 
illustrated  that  the  heat  transfer  characteristics  of  the  nanofluid 
can  be  enhanced  as  the  mean  nanoparticle  diameter  is  decreased 
from  250  to  5  nm.  They  attributed  these  phenomena  to  the 
dominant  effect  of  the  Brownian  motion  caused  by  heat  convec¬ 
tion.  However,  the  heat  transfer  performance  of  the  nanofluid  was 
found  to  be  less  significant  as  the  dimensionless  total  thermal 
conductivity  of  the  nanofluid  is  close  to  unity  due  to  the  increase 
of  nanoparticle  sizes.  This  contradictory  effect  of  nanofluids 
was  explained  to  be  mainly  caused  by  the  effective  dynamic 
viscosity. 

Abu-nada  and  Chamkha  [42]  focused  on  the  study  of  the 
natural  convection  heat  transfer  characteristics  in  a  differentially 
heated  enclosure  filled  with  Cuo-Eg-water  nanofluid.  Overall,  the 
study  indicated  that  the  effects  of  the  viscosity  models  are 
predicted  to  be  more  predominant  on  the  behavior  of  the  average 
Nusselt  number  than  the  influence  of  the  thermal  conductivity 
models.  Also,  the  enclosure  aspect  ratio  is  predicted  to  have 
significant  effects  on  the  behavior  of  the  average  Nusselt  number 
which  decreases  as  the  enclosure  aspect  ratio  increases. 

Abu-nada  [43]  presented  a  solution  of  Rayleigh-Benard  convec¬ 
tion  problem  of  CuO-water  nanofluids  flowing  in  a  rectangular 
enclosure.  The  results  were  presented  over  a  wide  range  of  Rayleigh 
numbers  and  volume  fractions,  using  variable  thermal  conductivity 
and  variable  viscosity  models.  It  has  been  shown  that  For  Ra  >  103, 
the  average  Nusselt  number  was  reduced  by  increasing  the  volume 
fraction  of  nanoparticles.  However,  for  Ra=103,  the  average  Nusselt 
number  was  enhanced  by  increasing  the  volume  fraction  of  nano¬ 
particles.  Abu-Nada  explained  such  deterioration  in  heat  transfer 
observed  experimentally  and  showed  that  the  influence  of  nano¬ 
particles  elucidates  two  opposing  effects  on  the  energy  transport:  a 
favorable  effect  driven  by  the  presence  of  high  thermal  conductive 
nanoparticles  and  undesirable  effect  promoted  by  high  level  of 
viscosity  experienced  at  high  volume  fraction  of  nanoparticles, 
where  at  high  Rayleigh  number  the  undesirable  effect  of  high 
viscosity  becomes  dominant,  which  causes  deterioration  in  heat 
transfer.  The  results  obtained  using  variable  properties  of  nanofluid 
were  compared  with  those  based  on  constant  property  simulations 
and  it  was  found  that  the  temperature  influence  is  small  compared 


to  the  influence  of  high  viscosity  brought  by  the  presence  of  high 
concentration  of  nanoparticles. 

Corcione  [44]  presented  two  general  correlations  for  thermal 
conductivity  and  viscosity  of  nanofluids  using  many  available 
experimental  data.  These  two  correlations  are  sensitive  to  types 
of  the  base  fluid  and  nanoparticle,  diameter  of  nanoparticle  and 
the  temperature  of  the  nanofluid.  Using  these  new  correlations, 
Corcione  [45]  investigated  the  heat  transfer  features  of  buoyancy 
driven  nanofluid  inside  rectangular  enclosures  differentially 
heated  at  the  vertical  walls.  They  assumed  that  nanofluids  behave 
more  like  single-phase  fluid  rather  than  a  conventional  solid- 
liquid  mixture,  thus,  they  extended  the  heat  transfer  correlations 
available  in  literature  for  single  phase  natural  convection  in 
enclosure  to  nanoparticle  suspension,  by  replacing  the  thermo¬ 
physical  properties  appearing  in  them  with  the  nanofluid  effec¬ 
tive  properties  calculated  at  the  average  temperature.  Results 
showed  the  existence  of  an  optimal  particle  loading  for  maximum 
heat  transfer.  Specifically,  for  any  assigned  combination  of  solid 
and  liquid  phases,  the  optimal  volume  fraction  was  found  to 
increase  slightly  with  decreasing  the  nanoparticle  size,  and  to 
increase  much  more  remarkably  with  increasing  both  the  nano¬ 
fluid  average  temperature  and  the  slenderness  of  the  enclosure. 

Ghasemi  and  Aminossadati  [46]  presented  the  results  of  a 
numerical  study  on  the  natural  convection  in  a  right  triangular 
enclosure,  with  a  heat  source  on  its  vertical  wall  and  filled  with  a 
water-Cuo  nanofluid.  The  viscosity  and  the  thermal  conductivity 
of  the  nanofluid  were  determined  based  on  the  models  presented 
by  Koo  and  Kleinstreuer  [47,48]  in  which  a  Brownian  motion  term 
was  added  and  valid  in  the  range  of  1  <  <fi  <  4%.  It  was  noticed 
that  when  Brownian  motion  is  neglected,  the  heat  transfer  rate 
continuously  increases  with  increasing  the  solid  volume  fraction 
at  all  Rayleigh  numbers,  whereas,  at  high  Rayleigh  numbers,  an 
optimum  solid  volume  fraction  can  be  found,  which  results  in  the 
maximum  heat  transfer  when  Brownian  motion  is  considered. 

Santra  et  al.  [49]  simulated  the  behavior  of  heat  transfer  due  to 
natural  convection  in  a  differentially  heated  cavity.  They  calcu¬ 
lated  the  shear  stresses  using  Ostwald-de  Waele  model  [50]  for 
an  incompressible  non-Newtonian  fluid.  It  was  observed  that  the 
heat  transfer  rate  decreases  with  increase  in  solid  volume  fraction 
for  a  particular  Rayleigh  number.  However,  it  increases  with 
Rayleigh  number  for  a  particular  solid  volume  fraction. 

Rashmi  et  al.  [51]  performed  a  numerical  study  of  natural 
convection  for  Al203-water  nanofluid  inside  a  horizontal  cylinder, 
using  FLUENT  commercial  code.  Stationary  particle  model  was 
employed  for  thermal  conductivity.  The  numerical  results  showed 
decrease  in  heat  transfer  with  increase  in  particle  volume  frac¬ 
tion.  Numerical  predictions  were  also  presented  to  be  compared 
with  experimental  results  of  Putra  et  al.  [76].  The  comparison 
showed  a  good  agreement  with  the  experimental  findings. 

Das  and  Ohal  [52]  conducted  a  numerical  investigation  of 
natural  convection  inside  a  square  cavity  containing  Cu-water 
nanofluid  with  half-active  and  half-insulated  vertical  walls.  They 
considered  five  different  cases  based  on  the  position  of  the  active 
walls.  It  was  found  that  the  nanoparticles  when  immersed  in  a 
fluid  are  capable  of  increasing  the  heat  transfer  capacity  of  the 
base  fluid  and  this  effect  is  more  pronounced  for  high  solid 
volume  fraction. 

Ghasemi  et  al.  [53]  investigated  the  effects  of  magnetic  field  on 
natural  convection  in  a  nanofluid-filled  square  enclosure  for 
different  Hartmann  number,  Rayleigh  number  and  nanoparticle 
volume  fraction.  They  indicated  that  the  heat  transfer  rate 
increases  with  an  increase  of  the  Rayleigh  number  but  it 
decreases  with  an  increase  of  the  Hartmann  number.  An  increase 
of  the  solid  volume  fraction  may  result  in  enhancement  or 
deterioration  of  the  heat  transfer  performance  depending  on  the 
value  of  Hartmann  number  and  Rayleigh  numbers. 
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Oztop  et  al.  [54]  studied  the  problem  of  steady  state  natural 
convection  in  an  enclosure  filled  with  different  types  of  nano¬ 
particles  (A1203  and  Ti02)  numerically.  It  was  observed  that  the 
addition  of  nanoparticles  into  water  affects  the  fluid  flow  and 
temperature  distribution  especially  for  higher  Rayleigh  numbers. 
In  addition,  an  enhancement  in  heat  transfer  rate  was  registered 
for  the  whole  range  of  Rayleigh  numbers.  However,  low  Rayleigh 
numbers  showed  more  heat  transfer  increment  compared  to  high 
Rayleigh  numbers.  It  was  also  noticed  that  the  difference  between 
A1203  and  Ti02  enhancement  becomes  more  pronounced  at  high 
volume  fractions  of  nanoparticles. 

Aminossadati  and  Ghasemi  [55]  performed  a  numerical  study 
of  natural  convection  in  a  square  cavity  filled  with  a  water-CuO 
nanofluid.  They  considered  two  pairs  of  heat  source-sink  to  cover 
the  entire  length  of  the  bottom  wall  of  the  cavity  while  the  other 
walls  were  thermally  insulated.  The  results  showed  that  regard¬ 
less  of  the  position  of  the  pairs  of  source-sink,  the  heat  transfer 
rate  increases  with  an  increase  of  the  Rayleigh  number  and  the 
solid  volume  fraction. 

Jahanshahi  et  al.  [56]  examined  the  influence  of  uncertainties 
due  to  adopting  various  formulas  for  the  effective  thermal  con¬ 
ductivity  of  Silica-water  on  heat  transfer  characteristics  for 
natural  convection  in  a  square  enclosure.  An  experimental  setup 
was  used  to  extract  the  conductivity  value  of  nanofluid.  They 
predicted  the  effective  thermal  conductivity  of  the  nanofluid  using 
the  model  proposed  by  Hamilton  et  al.  [57]  and  experimental 
results.  Their  findings  showed  a  considerable  enhancement  in  heat 
transfer  by  increasing  the  solid  volume  fraction  at  any  Rayleigh 
number  using  experimental  thermal  conductivity,  whereas  a 
decrease  in  heat  transfer  using  theoretical  thermal  conductivity. 

Kefayati  et  al.  [58]  reported  the  numerical  results  of  natural 
convection  in  enclosures  using  water/Si02  nanofluid.  The  main  aim 
of  their  study  was  to  identify  the  ability  of  Lattice  Boltzmann 
method  (LBM)  for  solving  nanofluid  and  various  geometries.  The 
simulations  were  performed  for  different  Rayleigh  numbers,  solid 
volume  fraction  and  aspect  ratio  of  enclosure.  Experimental  data 
from  Jahanshahi  et  al.  [56]  were  used  to  predict  thermal  conductiv¬ 
ity  of  Si02-water  nanofluid.  It  was  concluded  that  the  Nusselt 
number  increases  with  volume  fraction  for  the  whole  range  of 
Rayleigh  numbers  and  aspect  ratios.  Also  the  effect  of  nanoparticles 
on  heat  transfer  augments  as  the  enclosure  aspect  ratio  increases. 

Aminossadati  and  Ghasemi  [59]  presented  a  numerical  inves¬ 
tigation  to  study  the  natural  convection  in  an  isosceles  triangular 
enclosure  with  a  heat  source  located  at  its  bottom  wall  and  filled 
with  an  Ethylene  Glycol-Cu  nanofluid.  The  results  indicated  that 
the  thermal  performance  of  the  enclosure  is  improved  with  an 
increase  in  the  Rayleigh  number  and  solid  volume  fraction.  The 
results  also  indicated  that  the  variation  of  heat  transfer  rate  with 
respect  to  the  enclosure  apex  angle,  heat  source  position  and 
dimensions  is  different  at  low  and  high  Rayleigh  numbers.  More¬ 
over,  the  results  obtained  from  the  modified  and  original  Maxwell 
models  indicated  that  the  heat  transfer  rates  obtained  based  on 
the  modified  Maxwell  model  are  generally  higher  than  those 
obtained  based  on  the  original  Maxwell  model. 

Shahi  et  al.  [60]  analyzed  three-dimensional  flow  and  heat 
transfer  in  a  single  ended  tube  with  non-uniform  heat  input.  They 
assumed  that  the  sealed  end  of  tube  to  be  adiabatic  and  also  the 
tube  opening  to  be  subjected  to  copper-water  nanofluid.  The  results 
indicated  that  the  maximum  overall  mean  Nusselt  number  is 
obtained  at  <j)  =  35°,  while  the  maximum  output  mass  flow  rate  is 
increasing  function  of  the  inclination  angle.  It  was  also  indicated 
that  both  the  Nusselt  number  and  maximum  output  mass  flow  rate 
are  increasing  function  of  solid  concentration,  but  the  presence  of 
nanoparticles  is  more  effective  at  the  smaller  inclination  angles. 

Ghasemi  and  Aminossadati  [61]  conducted  a  numerical  study  to 
investigate  the  natural  convection  cooling  of  an  oscillating  heat 


source  embedded  on  the  left  wall  of  an  enclosure  filled  with  a 
water-based  nanofluids  (Cu,  A1203,  Ti02).  They  examined  the  effects 
of  Rayleigh  number,  solid  volume  fraction,  heat  source  position,  type 
of  nanofluids  and  oscillation  period  on  the  performance  of  the 
enclosure.  The  study  showed  that  the  addition  of  nanoparticles  into 
the  pure  water  improves  its  thermal  conductivity  and  enhances  the 
heat  removal  from  the  heat  source.  At  low  Rayleigh  number,  this 
effect  is  more  noticeable  for  copper  nanoparticles.  It  was  also 
showed  that  Cu  and  Ti02  are  the  most  and  the  least  effective 
nanoparticles,  respectively,  in  the  heat  removal  process. 

Elhajjar  et  al.  [62]  used  laws  of  thermodynamics  for  specific 
heat  capacity  and  thermal  expansion  coefficient  to  study  numeri¬ 
cally  the  Rayleigh-Benard  convection  for  three  types  of  nano¬ 
fluids  (water-based  A1203,  Cu  and  CuO).  It  was  observed  that  for  a 
fixed  value  of  the  fluid  Rayleigh  number,  the  nanofluid  Rayleigh 
number  decreases  with  the  volume  fraction  of  nanoparticles, 
meaning  that  the  nanoparticles  delay  the  onset  of  convection. 

Sheikhzadeh  et  al.  [63]  numerically  modeled  the  buoyancy- 
driven  fluid  flow  and  heat  transfer  in  a  square  cavity  with 
partially  active  side  walls  filled  with  Cu-water  nanofluid.  They 
studied  the  effects  of  the  position  of  active  portions  of  the  side 
walls,  and  volume  fraction.  The  result  showed  that  the  average 
Nusselt  number  increases  with  increasing  both  the  Rayleigh 
number  and  the  volume  fraction.  Moreover,  the  maximum  aver¬ 
age  Nusselt  number  for  the  high  and  the  low  Rayleigh  numbers 
occur  for  the  bottom-middle  and  the  middle-middle  locations  of 
the  thermally  active  parts,  respectively. 

Saleh  et  al.  [64]  studied  the  heat  transfer  enhancement 
utilizing  water-Cu  and  water-Al203  nanofluids  in  a  trapezoidal 
enclosure.  Various  inclination  angles  of  the  sloping  wall,  solid 
volume  fractions,  and  Grashof  numbers  were  considered  and  the 
flow  and  temperature  fields  as  well  as  the  heat  transfer  was 
analyzed.  They  found  that  the  structure  of  the  fluid  flow  within 
the  enclosure  depends  upon,  Grashof  number,  inclination  angle  of 
sloping  wall  and  nanoparticles  concentration  and  type.  In  addi¬ 
tion,  the  Cu  nanoparticles  with  high  volume  fraction  combines 
with  an  acute  sloping  wall  are  most  effective  in  enhancing 
performance  of  heat  transfer  rate. 

Mahmoodi  [65]  numerically  investigated  the  free  convection 
fluid  flow  and  heat  transfer  of  various  water  base  nanofluids.  The 
left  and  right  walls  of  the  cavity  are  maintained  at  constant 
temperature,  while  its  top  and  bottom  walls  are  insulated.  A  thin 
heater  is  located  inside  the  cavity  which  its  location  and  length. 
They  investigated  the  effects  of  the  Rayleigh  number,  position  and 
location  of  the  heater,  the  volume  fraction  of  the  nanoparticles, 
and  various  types  of  nanofluids  on  the  fluid  flow  and  heat 
transfer.  It  was  observed  that  for  all  cases,  the  rate  of  heat 
transfer  increases  with  increase  in  Rayleigh  number  and  the 
volume  fraction  of  the  nanoparticles.  It  was  also  observed  that 
at  low  Rayleigh  numbers  the  type  of  nanofluid  does  not  affect  the 
heat  transfer  rate  while  at  high  Rayleigh  numbers  maximum  rate 
of  heat  transfer  is  obtained  by  Ag-water  nanofluid  and  its 
minimum  is  obtained  by  Ti02-water  nanofluid. 

Lai  and  Yang  [66]  performed  a  numerical  study  to  examine  the 
natural  convective  thermal  efficacy  of  Al203/water  nanofluid  for 
two  circumstances,  i.e.  at  fixed  Rayleigh  number  and  at  fixed 
temperature  difference  across  the  enclosure,  using  the  lattice 
Boltzmann  method.  Results  indicated  that  the  average  Nusselt 
number  increases  with  the  increase  of  Rayleigh  number  and 
nanoparticle  volume  concentration.  In  addition,  the  relative 
Nusselt  number  increases  with  an  increasing  average  fluid  tem¬ 
perature  and  the  increase  of  temperature  difference  imposed 
between  the  side  walls.  The  deviations  between  the  predicted 
values  of  thermal  conductivity  using  various  models  were  found 
not  to  be  remarkable,  while  significant  under-estimations  of 
nanofluid  viscosity  were  revealed  using  different  models. 
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3.2.  Two  phase  model 

Two-phase  model  was  proposed  to  get  realistic  results.  Tzou 
[67,68]  analyzed  analytically  the  onset  of  convection  in  a  horizontal 
layer  uniformly  heated  from  below  (the  Rayleigh-Benard  problem). 
They  used  a  model  proposed  by  Buongiorno  [11]  for  convective 
transport  in  nanofluid  incorporating  the  effects  of  Brownian  diffusion 
and  thermophoresis.  In  the  study  of  Tzou  [67],  the  case  of  stress-free 
boundaries  was  analyzed.  In  his  other  work  [68],  the  combination  of 
one  free  and  one  rigid  boundary  was  investigated,  and  the  case  of  two 
rigid  boundaries  was  briefly  mentioned.  On  the  basis  of  his  calcula¬ 
tions,  Tzou  [67,68]  concluded  that  for  nanofluids  the  critical  Rayleigh 
number  (the  value  of  the  Rayleigh  number  at  which  instability 
appears)  was  lower  by  one  or  two  orders  of  magnitude  than  for 
regular  fluids. 

The  Cheng-Minkowycz  problem  was  solved  by  Nield  and 
Kuznetsov  [69]  for  natural  convective  boundary-layer  flow  in  a 
porous  medium  saturated  by  a  nanofluid.  They  included  Brow¬ 
nian  and  thermophoresis  effects  and  similarity  solution  was 
presented.  In  another  study  [70],  they  made  an  analytical  study 
at  the  onset  of  convection  in  a  horizontal  layer  of  a  porous 
medium  saturated  by  a  nanofluid  by  including  Brownian  and 
thermophoresis.  Again,  Nield  and  Kuznetsov  [71]  presented  a 
work  on  a  linear  stability  analysis  for  the  onset  of  natural 
convection  in  a  horizontal  nanofluid  layer.  The  similar  case  is 
studied  for  natural  convective  boundary-layer  flow  of  a  nanofluid 
past  a  vertical  plate  by  Kuznetsov  and  Nield  [72].  Khan  and  Pop 
[73]  published  a  paper  on  boundary-layer  flow  of  a  nanofluid  past 
a  stretching  sheet  as  a  first  paper  in  that  field.  Their  model  used 
for  the  nanofluid  incorporates  the  effects  of  Brownian  motion  and 
thermophoresis.  They  have  taken  into  account  the  Pr  number, 
Lewis  number,  Brownian  motion  numbers,  Nb  and  thermophor¬ 
esis  number,  Nt.  They  indicated  that  the  reduced  Nusselt  number 
is  a  decreasing  function  of  each  dimensionless  number  and 
Sherwood  number  is  an  increasing  function  of  higher  Prandtl 
number  and  a  decreasing  function  of  lower  Pr  number  for  each  Le, 
Nb  and  Nt  numbers.  Recently,  Pakravan  and  Yaghubi  [74]  studied 
the  thermophoresis,  Brownian  and  Dufour  effect  on  natural 
convective  heat  transfer  of  nanofluids  simultaneously.  They 
showed  that  the  effects  of  Dufour  strongly  decreases  with 
temperature  and  Nusselt  number  of  natural  convection  of  nano¬ 
fluids  increases  with  mean  temperature  of  mixture. 

Recently,  Haddad  et  al.  [75]  studied  numerically  the  natural 
convection  heat  transfer  and  fluid  flow  of  CuO-water  nanofluid  in 
an  open  cavity  heated  from  the  bottom.  A  two  component  non- 
homogenous  equilibrioum  model  was  used  for  the  nanofluid, 
incorporating  the  effects  of  Brownian  motion  and  thermophoresis. 
Variable  thermal  conductivity  and  variable  viscosity  were  taken 
into  account.  It  was  showed  that  higher  heat  transfer  is  formed  with 
the  presence  of  Brownian  and  thermophoresis  effect.  It  was  also 
found  that  heat  transfer  decreases  with  increasing  of  nanoparticle 
volume  fraction  due  to  increasing  of  viscosity  of  the  nanofluid. 

4.  Experimental  investigation 

Experimental  works  on  natural  convection  in  nanofluid  enclosure 
are  very  limited  due  to  difficulties  of  preparation  of  nanofluid  and 
measurements  of  the  parameters.  Putra  et  al.  [76]  made  an  experi¬ 
mental  work  on  natural  convection  of  A1203  and  CuO-water  nano¬ 
fluid  inside  a  horizontal  cylinder  heated  and  cooled  from  the  two 
ends  respectively.  They  found  that  the  presence  of  nanaoparticles 
deteriorates  the  natural  convective  heat  transfer  systematically.  They 
observed  a  systematic  degradation  of  natural  convective  heat  transfer 
with  increasing  particle  concentration  characterized  by  decreasing 
Nusselt  number  for  a  given  Rayleigh  number.  The  decrease  was 


observed  to  be  more  severe  for  aspect  ratio  of  1.0  than  that  of  0.5. 
They  ascribed  the  possible  reasons  of  deterioration  to  the  effects  of 
particle-fluid  slip  and  sedimentation  of  nanoparticles. 

Wen  and  Ding  [77]  conducted  an  experimental  study  of  natural 
convective  heat  transfer  of  Ti02  (30-40  nm)/water  nanofluid  in  a 
vessel  which  was  composed  of  two  horizontal  aluminum  disks  of 
diameter  240  nm  and  thickness  10  nm  separated  by  a  10  nm  gap. 
The  upper  surface  of  the  top  aluminum  disk  was  exposed  to  the 
open  space  and  cooled  by  compressed  air.  They  investigated  both 
the  transient  and  steady  heat  transfer  coefficient  for  different 
concentrations  of  nanofluids.  They  found  that  nanofluids  decrease 
the  natural  convective  heat  transfer  coefficient.  Furthermore,  such 
deterioration  increases  with  nanoparticle  concentrations.  They 
attributed  the  reason  of  their  observation  to  the  particle-particle 
interaction  and  modifications  of  the  dispersion  properties. 

Nnanna  [78]  presented  an  experimental  system  to  investigate 
the  natural  convective  heat  transfer  in  a  differentially  heated 
enclosure,  filled  with  Al203-water.  It  was  observed  for  small  volume 
fraction  0.2  <  </>  <  2%  the  presence  of  nanoparticles  does  not  impede 
the  free  convective  heat  transfer,  rather  it  augments  the  rate  of  heat 
transfer.  However,  for  large  volume  fraction  c/)  >  2%,  the  convective 
heat  transfer  coefficient  declines  due  to  reduction  in  the  Rayleigh 
number  caused  by  increase  in  kinematic  viscosity. 

Li  and  Peterson  [79]  performed  an  experimental  investigation  on 
the  natural  convection  heat  transfer  characteristics  of  Al203/water 
nanofluids  in  a  vertical  cylinder  enclosure  heated  from  below. 
Compared  with  the  base  fluid  (distilled  water),  the  heat  transfer 
rate  across  the  enclosure  was  found  increasingly  deteriorated  with 
the  volume  fraction  of  the  nanoparticles  in  the  nanofluid.  To  further 
explain  the  possible  causes,  other  than  the  relatively  enhanced 
viscosity  of  the  nanofluid,  for  the  deteriorated  heat  transfer  efficacy 
of  using  the  nanofluid,  a  flow-  visualization  experiment  of  sub¬ 
micron  polystyrene-water  suspension  in  rectangular  enclosure  was 
performed  to  infer  the  detrimental  influences  due  to  the  Brownian 
motions  and  thermophoresis  movements  of  the  particles. 

Ho  et  al.  [80]  reported  experimental  results  which  illustrated  the 
natural  convective  heat  transfer  of  Al203-water  nanofluid  in  a 
vertical  square  enclosure.  The  thermophysical  properties  including 
the  particle  size,  dynamic  viscosity,  density  and  thermal  conductivity 
were  measured  as  a  function  of  temperature  as  well  as  of  the 
volumetric  particle  fraction  in  the  nanofluid.  The  experimental 
results  illustrated  systematic  heat  transfer  degradation  for  the 
nanofluid  containing  nanoparticles  of  </)  >  2%  over  the  entire  range 
of  the  Rayleigh  number  considered.  However,  for  the  nanofluid 
containing  much  lower  particle  fraction  of  0.1%,  a  heat  transfer 
enhancement  of  around  1 8%  compared  with  that  of  water  was  found 
to  arise  in  the  largest  enclosure  at  sufficiently  high  Rayleigh  number. 

Table  3  shows  the  summary  of  published  experimental  and 
theoretical  investigations  of  the  natural  convective  heat  transfer 
performance  of  various  nanofluids. 

5.  Recent  Studies 

Recently,  several  studies  emerged  that  analyze  natural  con¬ 
vection  heat  transfer  enhancement  using  nanofluids.  Mahmoodi 
[90]  proposed  a  numerical  investigation  to  describe  the  natural 
convection  fluid  flow  and  heat  transfer  of  Cu-water  nanofluid 
inside  L-shaped  cavities.  The  effects  of  the  Rayleigh  number,  the 
aspect  ratio  and  the  volume  fraction  were  investigated.  The 
results  indicated  that  the  average  Nusselt  number,  for  all  ranges 
of  cavities  aspect  ratio,  increases  with  increase  in  the  Rayleigh 
number  and  the  solid  volume  fraction  of  the  nanofluid.  Also,  it 
was  indicated  that  the  rate  of  heat  transfer  was  higher  for  the 
L-shaped  cavities  with  lower  aspect  ratio.  Fattahi  [91]  investi¬ 
gated,  numerically,  natural  convection  flows  in  a  cavity  subject  to 
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Table  3 

Summary  of  experimental  and  theoretical  studies  on  natural  convection  of  nanofluids. 


Authors  Configuration 


Particle  Base  Experimental  /  Parameters  Observations 

material  fluid  Numerical 

method 


Khanafer  et  al. 
[17] 


Jang  and  Choi 

[21] 


Th 


\  \  \  \  \  \  \  \  \  \ 


H 


KKWWWW 

Cooled 


Kim  et  al.  [22] 


Jou  and  Tzeng 
[24] 


Rayfeigh-  Benard  colls  Temperature  profile 


Differentially  heated  cavity 


Polidori  et  al.  Heated  vertical  plate 
[25] 


Cu  Water 


Cu  Water 


Cu/  Water 

Diamond 


Cu/Ag  Water 


Cu  Water 


y-  A1203  Water 


Finite  volume 
method 


Finite-volume 

approach 


Theoretical  study 


Theoretical  study 


Finitel  difference 
method 


Theoretical  study 


1 03  <  qm  <  1 06 
0  <  (/>  <  5% 


103  <  Gr  <  105 


0  <  </>  <  25% 
dCu  =  10  nm 
A  =  1.0 

300  <T<  325  K 
dCu  =  6vnm 

rfdiamond  =  2  nm 

</>  <  5% 

0  <  (j)  <  20% 


1 03  <  Ra  <  1 06 

0  <  cp  <  20% 

dCu  =  10  nm 
A  =  0.5,  1.0,  2.0 
0  <  </>  <  4% 
dy-  A1203  =  42  nm 


Hwang  et  al.  [26] 


Ho  et  al.  [29] 


Oztop  and  Abu- 
nada  [30] 


Ogiit  [31] 


Abu-nada  and 
Oztop  [32] 


A1203  Water 


A1203  Water 


A1203/  Water 
Cu/Ti02 


Cu/  Water 

A1203/ 

CuO/Ag  / 

Ti02 


Cu  Water 


Theoretical  study 


Finite  volume 
method 


Finite  volume 
method 


PDQ.  method 


Finite  volume 
method 


300  <T<  325  K 
10  <  d  <  50  nm 

(/><  5% 

1 03  <  Ra  <  1 06 

0  <  4>  <  4% 

A  =  1.0 


1 03  <  Ra  <  5  x  105 
0  <  </>  <  20% 

0.5  <  A  <  2 


0  <  (j)  <  20% 


1 04  <  Ra  <  1 06 

0  <  (p  <  90° 

A  =  1.0 

0°<(p<  120° 

1 03  <  Ra  <  1 05 

0<</)<  10% 
0<  (p  <  120° 

A  =  1 


Both  Nusselt  number  and  maximum  output 
mass  flow  rate  are  increasing  function  of 
solid  concentration. 


The  nanofluid  heat  transfer  rate  increases 
with  an  increase  in  the  nanoparticles 
volume  fraction. 


Nu  =  0.5163(0.4436  +  (/)1  0809)Gr03123 


The  nanofluids  are  more  stable  and  have 
larger  heat  transfer  coefficient  compared 
with  a  classical  fluid. 


Heat  transfer  coefficient  of  a  nanofluid  is 
enhanced  by  all  parameters  with  respect  to 
the  volume  fraction. 


Increasing  the  volume  fraction  of  nanofluids 
cause  an  increase  in  the  average  heat 
transfer  coefficient. 

The  sensitivity  to  the  viscosity  model  used 
seems  undeniable  and  plays  a  key  role  in 
the  heat  transfer  behavior. 

The  heat  transfer  is  decreased  as  the  size  of 
nanoparticles  increases,  or  the  average 
temperature  of  nanofluids  is  decreased. 


Significant  difference  between 
enhancements  in  the  dynamic  viscosity 
estimated  from  the  two  adopted  formulas 
leads  to  contradictory  heat  transfer  efficacy 
of  the  nanofluid  Nu  =  C(  1  +cf))mRan 


The  heat  transfer  increases  with  increasing 
the  value  of  volume  fraction  of 
nanoparticles. 


The  average  heat  transfer  rate  increases 
significantly  as  particle  volume  fraction 
increases. 


The  effect  of  nanoparticle  concentration  on 
Nusselt  number  is  more  pronounced  at  high 
volume  fraction  than  at  low  volume 
fraction. 
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Table  3  ( continued ) 


Authors  Configuration 


Particle  Base  Experimental  /  Parameters  Observations 

material  fluid  Numerical 

method 


Abouali  and 
Falahatpisheh 
[33] 


Abu-nada  [34] 


AdlabMfc 


Aminossadati 
and  Ghasemi 


[37] 


Abu-nada  et  al. 
[38] 


Lin  and  violi  [39] 


Lrfl 

| 

(ftnodrl  ' 

j  i 


t 


Bccndin 
•  Ptnwk 


t 


He*n*r'Tg) 


2. 


Corcione  [45] 


« - w - - 

Differentially  heated  enclosure 


A1203  Water 


A1203  Water 


Cu/Ti02/  Water 
Al203/Ag 


A1203/  Water 
CuO 


A1203  Water 


CuO  Water/ 
EG 


CuO  Water 


Cu/Ti02/  Water/ 
A1203  EG 


Finite  volume 
method 


Finite  volume 
method 


Finite  volume 
approach 


Finite  volume 
method 


Weighting 
function  method 


Finite  volume 
method 


Finite  volume 
method 


Theoretical  study 


103  <  Gr  <  105 

0  <  (/>  <  6% 

1  <A<5 


1 03  <  Ra  <  1 05 

0  <  (f)  <  9% 

A  =  0.2,  0.4,  0.8 


1 03  <  Ra  <  1 0s 
0  <  (/>  <  20% 

A  =  1.0 


1 03  <  Ra  <  1 05 

0  <  (/>  <  9% 

1/2  <  A  <  2 
dc uo  =  29  nm 
dA 1203=47  nm 


103  <  Gr  <  106 
5  <  dAi203  <  250  nm 
0  <4>  <5% 

A  =  1.0 


1 03  <  Ra  <  1 05 

0  <  </>  <  6% 

1/2  <  A  <  2.0 


1 03  <  Ra  <  1 06 

0  <  0  <  9% 

A  =  2.0 


249<T<324I< 
25  <  d  <  150  nm 
1  <A<40 

<f><6% 


CuO 


Water  Finite  volume 
method 


1 03  <  Ra  <  1 0s 

\<(j)<4% 

1/2  <  A  <  2.0 


The  thermal  conductivity  of  Jang  and  Choi 
[22]  agrees  well  with  the  experiment  and 
shows  deterioration  of  natural  convection 
heat  transfer.  Nunf  =  NUf(  1  -ccf)) 


For  Ra  >  104,  the  average  Nusselt  number 
was  reduced  by  increasing  the  volume 
fraction  of  nanoparticles.  However,  for 
Ra  =  103,  the  average  Nusselt  number 
increased  by  increasing  the  volume  fraction 
of  nanoparticles. 


Adding  nanoparticles  into  pure  water 
improves  its  cooling  performance  especially 
at  low  Rayleigh  numbers. 


The  average  Nusselt  number  was  more 
sensitive  to  the  viscosity  models  than  to  the 
thermal  conductivity  models. 


Enhancement  and  mitigation  heat  transfer 
can  be  attributed  to  the  dominant  effect  of 
the  Brownian  motion  caused  by  heat 
convection. 


The  effects  of  the  viscosity  models  are  more 
predominant  on  the  behavior  of  the  average 
Nusselt  number  than  the  influence  of  the 
thermal  conductivity  models. 


For  Ra  >  103,  the  average  Nusselt  number 
was  reduced  by  increasing  the  volume 
fraction  of  nanoparticles.  However,  for 
Ra  =  103,  the  average  Nusselt  number  was 
enhanced  by  increasing  the  volume  fraction 
of  nanoparticles 

The  optimal  volume  fraction  is  found  to 
increase  slightly  with  decreasing  the 
nanoparticle  size,  and  to  increase  much 
more  remarkably  with  increasing  both  the 
nanofluid  average  temperature  and  the 
slenderness  of  the  enclosure. 

When  Brownian  motion  is  neglected,  the 
heat  transfer  rate  continuously  increases 
with  increasing  the  solid  volume  fraction  at 
ah  Rayleigh  numbers. 
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Table  3  ( continued ) 


Authors  Configuration 


Particle  Base  Experimental  /  Parameters  Observations 

material  fluid  Numerical 

method 


Santra  et  al.  [49] 


Rashmi  et  al.  [51  ] 


Das  and  Ohal 
[52] 


\\\\\\\\\\ 


adiabatic 


Tc 


Ti 


Cu  Water 


A1203  Water 


Cu  Water 


Finite  volume 
method 


Fluent  software 


Finite  volume 
method 


1 04  <  Ra  <  1 0 
0  <4>  <5% 
dCu  =  100  nm 
A  =  1.0 


cf)  =  l%,4% 

A  =  1.0 


104  <  Gr  <  107 
0  <  4>  <  20% 

A  =  1.0 


adiabatic 


Ghasemi  et  al. 
[53] 


Oztop  et  al.  [54] 


Aminossadati 
and  Ghasemi 


[55] 


Jahanshahi  et  al. 
[56] 


Kefayati  et  al. 
[58] 


Adiabatic 

Water+AIOj 

Nanofluid 

Tk 

1  T‘ 

'1 

J 

X 

Adiabatic 

L 


A1203  Water  Finite  volume 

method 


Cu/Ti02/  Water  Finite  volume 
A1203  method 


CuO  Water  Finite  volume 

method 


Si02  Water  Finite  volume 

method/ 

Experimental 


Si02 


Water  Lattice 

Boltzmann 

method 


1 03  <  Ra  <  1 07 
1 03  <  Ha  <  1 07 

0  <  (f)  <  6% 

A  =  1.0 


1 03  <  Ra  <  1 05 

0<(f)<  10% 

0  <  cp  <90  ° 

A  =  1.0 


103  <  Gr  <  106 

0  <(f)  <4% 

A  =  1.0 


1 05  <  Ra  <  1 07 

0  <  4>  <  4% 

A  =  1.0 

dsi02  =  12  nm 


1 03  <  Ra  <  1 05 

0  <  cj)  <  4% 

1/2  <  A  <  2.0 


The  heat  transfer  rate  decreases  with 
increase  in  solid  volume  fraction  for  a 
particular  Rayleigh  number.  However,  it 
increases  with  Rayleigh  number  for  a 
particular  solid  volume  fraction. 


The  rate  of  heat  transfer  decreases  with 
increase  in  particle  volume  fraction. 


The  rate  of  heat  transfer  increases  with  an 
increase  in  the  nanoparticles  volume 
fraction. 


The  heat  transfer  rate  increases  with  an 
increase  of  the  Rayleigh  number  but  it 
decreases  with  an  increase  of  the  Hartmann 
number. 


An  enhancement  in  heat  transfer  rate  was 
registered  for  the  whole  range  of  Rayleigh 
numbers. 


The  heat  transfer  rate  increases  with  an 
increase  of  the  Rayleigh  number  and  the 
solid  volume  fraction. 


The  mean  Nusselt  number  increases  with 
volume  fraction  for  the  whole  range  of 
Rayleigh  numbers.  Although  by  using  the 
theoretical  formulations  for  conductivity  no 
enhancement  was  observed. 


The  Nusselt  number  increases  with  volume 
fraction  for  the  whole  range  of  Rayleigh 
numbers  and  aspect  ratios. 
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Table  3  ( continued ) 


Authors  Configuration 


Particle  Base  Experimental  /  Parameters  Observations 

material  fluid  Numerical 

method 


Aminossadati 
and  Ghasemi 
[59] 


Ghasemi  and 
Aminossadati 
[61] 


Cu 


EG 


Finite  volume 
method 


Cu/Ti02/ 

AI2O3 


Water 


Finite  volume 
method 


Elhajjar  et  al.  [62]  Rayleigh-Benard  configuration  Cu/Ti02/  Water 


Sheikhzadeh 
et  al.  [63] 


Saleh  et  al.  [64] 


Mahmoodi  [65] 


Z.  Haddad  et  al. 
[75] 


Putra  et  al.  [76] 


Cu/Ti02/ 

Ag 


Water 


Finite  volume 
method 


T,  I 


r*  r. 


Lai  and  Yang  [66] 


Al203 


Water 


Lattice 

Boltzmann 

method 


CuO 


Water  Finite  volume 
method 


<  tank 


CuO  / 

ai2o3 


Water  Experimental 
study 


. — 

H 

1 03  <  Ra  <  1 0s 

0<4>  <5% 


1 03  <  Ra  <  1 0s 
0  <  4>  <  20% 

A  =  1.0 


0  <  Ra  <  1 04 
0  <4>  <8% 

1 03  <  Ra  <  1 0s 


0<  <()<  15% 


103  <  Gr  <  105 


0  <  4>  <  26% 

0°  <cp  <  26° 


1 03  <  Ra  <  1 0s 

0<4><9  % 

A  =  1.0 


1 03  <  Ra  <  1 06 

0  <  (f>  <  4% 

0.006  <  A T<  10.5k 
A  =  1.0 


1 03  <  Ra  <  1 0s 

0  <  (j>  <  9% 

A  =  2.0 


1 07  <  Ra  <  1 09 

cf)  =  1%,4% 
c/ai203  =  131.2  nm 
dCu0  =  87.3  nm 
A  =  0.5,  1.0 


1.  CyltndrKal  block 

2.  End  cover  hcaung  «urf»c< 

3.  F.nd  co»w  as  cooling  surface 

4  Cap 

S.  Resistance  heating  elements 


6  The  piston  shaft 
7.  Cooling  water  inlet  and  outlet 
8  Narrow  tube 
9.  Thermocouples 


The  thermal  performance  of  the  enclosure  is 
improved  with  an  increase  in  the  Rayleigh 
number  and  solid  volume  fraction. 


The  addition  of  nanoparticles  into  the  pure 
water  improves  its  thermal  conductivity 
and  enhances  the  heat  removal  from  the 
heat  source. 


The  nanofluid  Rayleigh  number  decreases 
with  the  volume  fraction  of  nanoparticles. 
The  average  Nusselt  number  increases  with 
increasing  both  the  Rayleigh  number  and 
the  volume  fraction. 

Nuavg  =  (Q.4778^1  0  +  0.3554)Ra025 


The  Cu  nanoparticles  with  high  volume 
fraction  combines  with  an  acute  sloping 
wall  are  most  effective  in  enhancing 
performance  of  heat  transfer  rate. 

/.  \  -11.707 

NuaVg  =  0.231GrjJ342/l'1 ' 


y<bfM 


Acute  sloping  wall  and  Cu  nanoparticles 
with  high  concentration  are  effective  to 
enhance  the  rate  of  heat  transfer. 


The  average  Nusselt  number  increases  with 
the  increase  of  Rayleigh  number  Ra  and 
nanoparticle  volume  concentration. 


Higher  heat  transfer  is  formed  with  the 
presence  of  Brownian  and  thermophoresis 
effect. 


A  systematic  degradation  of  natural 
convective  heat  transfer  with  increasing 
particle  concentration  for  a  given  Rayleigh 
number. 
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Table  3  ( continued ) 


Authors  Configuration 


Particle  Base  Experimental  /  Parameters  Observations 

material  fluid  Numerical 

method 


Wen  and  Ding 
[77] 


Nnanna [78] 


Experimental  setup:  (A)  filling  vessel.  (B>  discharging  vessel. 
(C)  top  aluminium  plate,  (D)  bottom  aluminium  plate.  (E) 
insulation  material.  (F)  insulation  cover.  (G)  level  indicator, 
(H)  heating  element.  (I)  heat  flux  sensor.  (J)  DC  power 
supply.  (K)  data  acquisition  system  (DAQh  (L)  computer. 


Ti02 


AI2O3 


Water  Experimental 
study 


Water  Experimental 
study 


dsi02  =  30-40  nm 
A  =  0.5,  1.0 
0<</><  1% 

1 04  <  Ra  <  1 06 


1  x  107<Ra<3  x  107 


0  <  </>  <  8% 
dA 1203  =  27  nm 

A  =  1.0 


Nanofluids  decrease  the  natural  convective 
heat  transfer  coefficient,  such  deterioration 
increases  with  nanoparticle  concentrations. 


Nusselt  increased  with  small  volume 
fraction, 0.2  <  </>  <  2%,  and  no  significant 
change  in  Nusselt  number  is  observed  in  the 
range  2%  <4>  <  7.9%. 

Nu  =  re-Ra^{e^ 


Li  and  Peterson 
[79] 


Ho  et  al.  [80] 


Mahmoodi  [90] 


Fattahi  [91] 


Kefayati  et  al. 
[92] 


Abouali  and 
Ahmadi  [93] 

Guo  et  al.  [94] 


Experimental  setup.) )  insulation.  2)  resistive  electrical  healer. 
3)  heated  wall.  4)  nanofluid  medium.  S)  thermocouple  wire. 
6)  the  black  dots  indicate  the  thermocouple  locations.  7) 
unhealed  wall.  8)  spout:  i 
=35  mm.  H= 215  mm. 


AI2O3 


Water  Experimental 
study 


F.xpcrimcntal  test  facility:  (l  eft)  I:  horizontal  plane  table:  2: 
base  plate:  3:  vertical  support:  4:  load  screw:  5:  heal 
dissipater:  6:  Thermal  insulation  wrapper:  7:  upper  copper 
bar;  8:  sample  charge  tubes;  9:  o-ring:  10:  lower  copper  bar: 
II:  heat  sink;  12:  coolant  inlet  and  outlet;  (Right)  the 
thermocouple  distribution  on  both  copper 
bars. 


A1203  Water  Experimental 

study 


V 


Differentially  heated  cavity 


Cu  Water  Finite  volume 

method  Simpler 
algorithm 


Cu/A1203  Water 


Lattice 

Boltzmann 

method 


Th 


Six  different  enclosure  shapes 


Cu 


CuO  / 
AI2O3 


Water  Lattice 

Boltzmann 

method 


Water  Fluent  software 


Differentially  heated  cavity 


AI2O3 


Water  Lattice 

Boltzmann 

method 


0.5%  <(f)<  6% 

dA 1203=47  nm 
8  x  103  <  Ra  <  2.8  x  104 


The  heat  transfer  rate  was  found 
increasingly  deteriorated  with  the  volume 
fraction  of  the  nanoparticles. 


6.21  x  105  <  Ra  <  2.56  x  108  Systematic  heat  transfer  degradation  for  the 
0.1%  <  4>  <  4%  nanofluid  containing  nanoparticles  of 

<f>  >  2%  over  the  entire  range  of  Rayleigh 
number. 


1 03  <  Ra  <  1 06 
0<</><  10  % 

A  =  0.2,0.4,0.6 


The  average  Nusselt  number  increases  with 
increase  in  the  solid  volume  fraction  of  the 
nanofluid. 


1 03  <  Ra  <  1 06 

0  <  </>  <  5% 


1 04  <  Ra  <  1 06 
0  <<f>  <5% 

1/2  <  A  <  2.0 


The  average  Nusselt  number  increases  for 
both  nanofluids.  Further,  the  effects  of  solid 
volume  fraction  for  Cu  are  stronger  than 
AI2O3. 

The  average  Nusselt  number  increases  with 
augmentation  of  Rayleigh  number  and  the 
volume  fraction  of  nanoparticles  for  whole 
ranges  of  aspect  ratios. 


1 03  <  Ra  <  1 05 
0<</><  15% 

A  =  1.0 


The  numerical  simulations  for  natural 
convection  of  nanofluids  could  be  obtained 
theoretically  in  a  much  simpler  way. 

The  critical  Rayleigh  number  is  lower  for 
nanofluid  than  that  for  base  fluid. 

Nu  =  0.05Ra°-45 
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Table  3  ( continued ) 


Authors  Configuration 


Particle  Base  Experimental  /  Parameters  Observations 

material  fluid  Numerical 

method 


Mahmoodi  and 
Hashemi  [95] 


.  * 


Cu 


Cu 


Water  Finite  volume  103<Ra<106 

method  Simpler  0<(/><10% 
algorithm  71  =  0.2,0.4,0.6,0.8 


Water  Finite  Element 
Method 


1 03  <  Ra  <  1 05 

0<4><6% 
45°<(p<  180  0 

A  =  1.0 


The  mean  Nusselt  number  increased  with 
increase  in  volume  fraction  regardless 
aspect  ratio  of  the  enclosure. 


The  effect  of  nanoparticles  is  more 
pronounced  at  low  Rayleigh  number  than  at 
high  Rayleigh  number. 


different  side  wall  temperatures.  The  fluid  in  the  cavity  was  a 
water-based  nanofluid  containing  A1203  and  Cu  nanoparticles. 
The  effective  thermal  conductivity  and  viscosity  of  nanofluid  were 
calculated  using  Chon  and  Brinkman  models.  The  results  indi¬ 
cated  that  by  increasing  solid  volume  fraction,  the  average 
Nusselt  number  increased  for  both  nanofluids.  Kefayati  et  al. 
[92]  numerically  studied  natural  convection  of  Cu-water  nano¬ 
fluid  in  an  open  enclosure  using  Lattice  Boltzmann  method  (LBM). 
The  results  indicated  that  the  average  Nusselt  number  increased 
with  augmentation  of  the  Rayleigh  number  and  the  volume 
fraction  of  nanoparticles  for  the  whole  range  of  aspect  ratios. 
However,  the  average  Nusselt  number  decreased  as  the  aspect 
ratio  increased.  Abouali  and  Ahmadi  [93]  analyzed  natural  con¬ 
vection  of  nanofluids  in  horizontal  and  titled  square,  horizontal 
and  vertical  annuli,  a  triangular  enclosure  and  the  Rayleigh- 
Benard  convection  configuration.  The  study  was  performed  for 
A1203  and  CuO  nanofluids.  The  results  were  compared  with  the 
prediction  of  a  theoretical  approach  which  uses  the  Nusselt 
correlations  available  in  literature  for  pure  fluids.  Guo  et  al.  [94] 
conducted  a  numerical  study  for  natural  convection  heat  transfer 
in  a  square  cavity  using  the  LBM.  They  derived  a  correlation  for 
the  mean  Nusselt  number.  Mahmoodi  and  Hashemi  [95]  pre¬ 
sented  a  numerical  study  of  natural  convection  fluid  flow  and 
heat  transfer  inside  C-shaped  enclosures  filled  with  Cu-water 
nanofluid  using  finite  volume  method.  It  was  found  that  the  mean 
Nusselt  number  increases  with  increase  in  the  Rayleigh  number 
and  volume  fraction  of  Cu  nanoparticles  regardless  of  aspect  ratio 
of  the  enclosure.  Also,  the  effect  of  nanoparticles  on  enhancement 
of  heat  transfer  at  low  Rayleigh  numbers  is  more  significant  than 
that  at  high  Rayleigh  numbers.  Moreover,  it  was  shown  that  the 
rate  of  heat  transfer  increases  as  the  C-shaped  enclosure  becomes 
narrower.  Soleimani  et  al.  [96]  proposed  the  control  volume  based 
finite  element  method  (CVFEM)  to  investigate  natural  convection 
heat  transfer  in  a  semi-annulus  enclosure  filled  with  Cu-water 
nanofluid.  The  inner  and  outer  semi  circular  walls  are  maintained 
at  constant  temperatures  while  the  two  other  walls  are  thermally 
insulated.  The  numerical  investigation  is  carried  out  for  different 
governing  parameters:  namely  the  Rayleigh  number,  nanoparticle 
volume  fraction  and  the  angle  of  turn  for  the  enclosure.  The 
effective  thermal  conductivity  and  viscosity  of  nanofluid  were 
calculated  using  the  Maxwell-Garnetts  (MG)  and  Brinkman  mod¬ 
els,  respectively.  It  was  noticed  that  the  streamlines  and  iso¬ 
therms  strongly  depend  on  the  volume  fraction  of  nanoparticle. 


6.  Conclusion 

This  paper  presented  a  comprehensive  review  for  the  recent 
work  published  on  heat  transfer  enhancement  in  natural  convec¬ 
tion  using  nanofluids.  The  paper  presented  experimental,  numer¬ 
ical  and  analytical  published  work  in  literature.  The  numerical 
work  included  single-phase  and  two-phase  models.  Various 
thermal  conductivity  and  viscosity  models  employed  by  various 
researchers  were  presented. 

The  literature  survey  shows  that  natural  convection  heat 
transfer  enhancement  using  nanofluid  is  still  controversial  and 
there  is  an  on-going  debate  on  the  role  of  nanoparticles  on  the 
heat  transfer  enhancement  since  there  is  a  lack  of  extensive 
studies  on  this  topic.  Clearly,  most  of  numerical  results  showed 
that  nanofluids  significantly  improve  the  heat  transfer  capability 
of  conventional  heat  transfer  fluid.  Whereas,  experimental  results 
showed  that  presence  of  nanoparticles  deteriorates  heat  transfer 
systematically.  It  was  observed  that  the  experimental  results 
were  mostly  obtained  using  one  type  of  nanofluid  (Al203-water). 
Thus,  benchmark  experiments  are  highly  desired  to  gage  the 
validity  of  the  numerical  results.  Concerning  the  numerical 
studies  which  try  to  explain  the  observed  anomalous  enhance¬ 
ment  in  heat  transfer,  it  will  be  necessary  to  consider  not  only  one 
single  phase  model  but  two-phase  approach,  which  seems  a 
better  model  to  describe  the  nanofluid  flow  since  slip  velocity 
between  the  particle  and  base  fluid  plays  important  role  on  the 
heat  transfer  performance  of  nanofluids. 
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